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ABSTRACT 


. Thinning of the Earth's lithosphere by heat advected to Its 
base Is a possible oechanlsn for continental rifting and continental and 
oceanic m1d*p1ate volcanism. It night also account for continental 
r1ft1ng>11ke processes and volcanism on Venus and Mars. Earth's 
continental lithosphere can be thinned to the crust In a few tens of million 
years by heat advected at a rate of 5 to 10 times the normal basal heat flux. 

This BHich heat Is easily carried to the lithosphere by mantle plumes. The 
continent is not required to rest over the mantle hot spot but may move at 
tens of mill loiters per year. Because of the constant level of crustal 
radioactive heat production, the ratio of the final to the Initial surface 
heat flow increases much less than the ratio of the final to initial basal 
l«at flow. For large increases in asthenospheric heat flow, the littosphere 
is almost thinned to the crust before any significant change in surface heat 
flow occurs. Uplift due to thermal expansion upon thinning is a few kilometers. 

The oceanic lithosphere can be thinned to the crust in less than 10 million years If 
the heat advection is at a rate around 5 or more times the basal heat flow into 
100 Ma old lithosphere. Uplift upon thinmng can compensate the subsidence of 
spreading and cooling lithosphere. The swells will rise to a common depth 
of about 4 km below sea-level if the uppermost mantle beneath the mid-oceanic 
ridge is partially melted to a slightly larger degree than the uppermost mantle 
beneath the swells. Uplift upon thinning on Mars is up to 10 km or more, thus 
explaining the height of volcanic regions like Tharsis. On Venus, uplift 
upon thinning would range from 0.5 - 4 km, in agreement with observed topographic 
variations. 
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INTRODUCTION 


On Earth, almost all volcanic activity Is confined to the plate margins. 
However, scattered around the globe on oceanic and continental plates are 
regions of Isolated volcanic activity called hot spots. These surface hot 
spots are manifestations of mantle hot spots. Insofar as geologically recent 
plate motions are concerned the hot spots are quasi -stationary with respect 
to the lower mantle (Minster et al., 1974; Solomon and Sleep, 1974; Xaula, 1975). 
Consequently, when lithospheric plates pass over mantle hot spots,* swells and 
sometimes trails of volcanoes are left behind. Examples on oceanic plates 
Include the Hawaiian and Emperor swells and Island and seamount chains, the 
Austral and Tuamotu swells and Islands, and the Iceland transverse ridge, (hi 
continents, hot spot tracks are more difficult to Identify because continental 
plates move slower and are thicker than oceanic plates. Morgan (1981) succeeded 
In extending hot spot tracks from the Atlantic into the An^rlcas and Africa. 

Hot spots on continents are often accompanied by dome or plateau uplift and 
continental rifting (see the hot spot map recently provided by Vogt (1981)). 
Prominent examples Include the East African dtmie and rift zone, the Tibesti 
Plateau, the German Rhelnische Schild and the Rheingraben Rift, the Rio Grande 
Rift and the Colorado Plateau, the Siberian Platform and the Baikal Rift. 

Tectonic structures that resemble continental domes, plateaus, and 
rifts on Earth are also found on Venus and Mars. On Venus, doming and rifting 
Is thought to have caused the uplifted Beta and Phoebe Reglones which have been 
compared to the topographically similar Kenyan and Ethiopian domes (Mai In and 
Saunders, 1977; Wood and Head, 1978; Masursky et al., 1980; Phillips et al., 

1981; McGill et al., 1981). Two more rift-type 1 Inears have been recognized 
from Ploneer-Venus radar altimetry by Schaber (1982); the "Aphrodite-Beta 
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linear" extending northeast .from Aphrodite Terra to Beta Reglo and the "Themis- 
Atia linear" trending northwest frcwi east of Themis Reglo to the northeastern 
tip of Aphrodite Terra (the location of Atla Reglo). Schaber (1982) suggests 
that Beta and Atla Reglones are dynamically supported volcanic terranes located 
at the Intersections of the currently active major Venusian rift systems. 
Vigorous volcanic activity has shaped the surface of Mars throughout most of 
Its history (Greely and Spudis, 1981). There a»*e large uplifted plateaus and 
huge shield volcanoes. The Tharsis plateau Is the site of Olympus Mons, perhaps 
the largest shield volcano In the solar systmn. There Is also a huge canyon 
complex. Valles Marinerls, that res«nbles the East African Rift system. 

It has been proposed that hot spots are the surface manifestations of 
mantle plumes, possibly of deep origin (Wilson, 1963a, b; Morgan, 1971). Other 
causes for hot spots may be thermal runaway In asthenospheric shear flow 
(Melosh and Ebel, 1979) or the subduction cf active ridges under continents 
(Damon, 1979). Yuen and Peltier (1980) propose that mantle plumes could rise 
from the unstable core-mantle boundary layer of a whole mantle convection cell. 
This origin of the plumes may also explain the geochemical peculiarities of hot 
spot volcanism (Hofmann and White, 1982; Spohn and Schubert, 1982a). According 
to Anderson (1981), plumes could originate in the upper 220 km of the mantle. 

He suggests that thermal blanketing by the thick continental lithosphere could 
heat and partially melt the upper asthenosphere causing upwelling of the 
partially melted rock. 

Steady-state models of mantle plumes have been developed by Parmentler 
et al. (1975) and Yuen and Schubert (1976). Plumes can have upward velocities 
as large as meters per year and advect heat at significant rates, up to 50 
times the basal heat flux into the lithosphere. Where a plume Impinges on the 
lithosphere the excess heat causes the lithosphere to thin. Heating and 
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replacement of lithosphere by hot plume material cause uplift due to thermal 
expansion and Isostatic equilibration. On continents, a dome may form and 
rifts may develop In the dome (Burke and Wilson, 1976; Neugebauer, 1978). 
Eventually, the continent may split along a line of connecting rift arms of 
several domes. It has been shown by Spohn and Schubert (1982b) that the 
continental lithosphere can be thinned to the crust on the time scale of tens 
of million years which Is usually associated with continental rifting. Other 
authors have also attributed continental rifting to lithospheric thinning by 
enhanced mantle heat flow (e.g. Crough and Thompson, 1976a; Gass et a1., 1978). 

As for the oceanic lithosphere. Detrick and Crough (1978) and Crough 
(1978) have shown that the subsidence histories of oceanic mid-plate swells 
can be explained by a model In which the lithosphere Is reheated and thinned 
as It passes over a mantle hot spot. However, using the thinning nwdel o^ 

Crough and Thompson (1976a) and a heat flux anomaly of 2 times normal, they 
found that It would take around 100 Ma to conduct enough heat Into the lithosphere 
to lift the Hawaiian swell to Its maximum observed height while the sv«ll 
originates In only a few million years. Our results are In accord with this 
conclusion (Fig. 5). Thinning rates from Crough and Thompson's (1976a) nwdel 
are typically even higher than those computed from our model (Spohn and 
Schubert, 1982b). Thus, higher heat flux, anomalies are required. We will shw 
that a factor of 10 Increase will suffice. Observations of conductive heat 
flow at the surface of the oceanic crust are not a good guide to the variations 
of heat flux at lithospheric depths. High heat flux anomalies at depth may 
never be observed as conductive heat flux at the surface. For one thing, it 
will take 20 Ha or more before the thermal anomaly reaches the surface because 
of conductive heat transfer across the thinning lithosphere. Furthermore, 
a considerable part of the advected energy from the plume will be consumed by 
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depth. A heat flux of 1 H m can be consumed by a layer of magma growing at 
the rate of 20 ma a~^. This heat will eventually reach the surface via vol- 
canlsm. Another efficient mode of heat extraction Is hydrothermal convect1(Wi 
from magmatic Intrusions. 

Turning to Mars and Venus, we note that uplift of the Tharsis region 
on Kars can be Isostatically compensated if tne crust Is underlain by a mantle 
of anomalous low density (Sleep and Phillips, 1979). Isostatic ccm^nsatlwi 
of uplift on Venus Is possible If asthenospheric density anomalies are kept 
fron relaxing (Phillips et al., 1^1). 

In this paper we apply our convective thinning model (Spohn and Schubert, 
1982b) to the oceanic and continental lithospheres of Earth and to the litho- 
spheres of Venus and Mars. We begin by briefly reviewing the model and Its 
general consequences. Some new results are presented for the continental 
lithosphere. We then proceed to discuss the convective thinning of the oceanic 
lithosphere and show that thinning to the crust can occur in only a few million 
years. We give an explanation of why swells rise to a depth of about 1800 m below 
the mid-oceanic ridge crest. Finally, we apply our model to Mars and Venus. 


THE MODEL 

The model consists of a horizontally infinite lithosphere of thickness 
1^, uniform density p, specific heat c, and thermal conductivity k. The lower 
boundary of the lithosphere, the lithosphere-asthenosphere boundary LAB, Is a 
rheological boundary located at a depth where the viscosity is 0(10 Pa s), 
the characteristic magnitude of the sublithospheric mantle viscosity (Cathles, 
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1975; Peltier, 1981). The viscosity of mofitle rock 1$ strorqly teiRpereture 
dependent and the LAB is conveniently marked by an isotherm or a given fraction 
of the mantle solidus. We assume that the LAB temperature varies linearly 
with depth z according to 

T^(z) • fR^z + . (1) 

It is also possible that the solidus is reached at the LAB. There is evidence 
for Incipient melting in the upper mantle (Anderson and Sammis, 1970). 

The upper layer of the lithosphere is the crust of thickness z^. Rad- 
iogenic heat production is assumed to be confined to this layer; it varies 
with depth according to 


A(z) = Ajj exp(- 0 -) , (2) 


where A is the surface volumetric heat production rate and 0 is a character- 
istic depth for the near surface concentration of radioactive elements (Lach- 
enbruch and Sass, 1977). initially, the lithosphere is at thermal equilibrium 
with the temperature T^(z) given by 

TqCz) • X * ^ T . 0 < 2 < 

' T • *0 - * - *0 • • 

where is the basal heat flux from the asthenosphere. The quantities q., 

A^, and D are related to the surface heat flux q^^ by 

**si “ **i ^0 * 

It can be easily demonstrated from ( 1 ) and ( 2 ) that 
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Is the temperature rise across the crust due to radiogenic heat production. 

At t « 0 we suddenly Increase the heat flux from the asthenosphere 
from to q by placing a plume beneath It, for example. Thermal equilibrium 
Is perturbed and the LAB starts to rise to a new equilibrium position corre- 
sponding to a final lithospheric thickness given by 



(7) 


The space left by the boundary will be Immediately filled by the plume. The 
rising lithosphere Is underlain by a stagnation point boundary layer (Schllchting, 
I960) across which temperature increases by Tp - T^. It Is also assumed that 
Tp Is a linear function of depth 


T_ = m_2 + T„. 
p p ps 


( 8 ) 


We study the time-dependent position of the LAB along the plume center- 

3 T 

line X = 0, where x is horizontal distance. Here, ~ 0 both within the plume 

2 7 

3^T 3^T 

and the stagnation point boundary layer. Also, — — y within the lithosphere, 

3x^ 3z^ 

Thus, the LAB at x = 0 will be horizontal and its upward velocity will not depend 
on derivatives with respect to x. Heat is transported across the thinning 
lithosphere purely by conduction. Temperature along x * 0 in the lithosphere 
therefore satisfies 


pc|I = k ^ + A(z) 
3z^ 


( 9 ) 


and is subject to the initial and boundary conditions 
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where 2^(t) 1$ the position of the LAB at time t. 

Energy conservation at the LAB provides an equation for Its inotlon. 

If the base of the lithosphere moves upward a distance 6z^ In time 6t, then the 
quantity of heat - pc(Tp(Zj^) - pe«f* unit horizontal area must be 

provided to allow the plume with temperature Tp to follow the uoward motion 
of the lithosphere base. If the LAB Is at the solidus and there Is some 
partial melt, an amount of heat > pL£z^ per unit horizontal area must also 
be supplied, where L Is the latent heat. This energy Is provided by the 
difference between the upward heat flux Into the base of the lithosphere q 
and the upward heat flux out of the base of the lithosphere . Thus 


we can write 


L ^ c(Tp(z^) 


T.(z,)) 


■3T 





. ( 11 ) 


Equations (9) through (11) constitute a nonlinear boundary value problem which 
we solve numerically using a finite difference scheme. Convergence has been 
tested by decreasing the spacing of the grid points. The Courant criterion 
has been used to guarantee stability of the solutions. 

Uplift due to thermal expansion of the thirning lithosphere and replace- 
ment of lithosphere by hot plume material can be calculated fran cur model 
by balancing the masses In lithospheric columns before and after thinning. After 
the LAB has reached its final equilibrium position Isostatic uplift in continen- 
tal regions Is given by 


Where is the expansion of the thinned lithosphere due to heating 



j (^ “ <1^; 

1T~ 


f 


(13) 
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and o Is the volumetric coefficient of thermal expansion. Themrerage density 
difference App between lithosphere prior to heating and plume material In 
lo - If Is given by 


where Is the average density of the Initial lithosphere, Is the 
density change due to partial melting, and 

t* s (t^ + 1^) . (15), 

The Isostatic uplift In oceanic regions Is 


u 


PQ * % 
f>o * % * 


Pw 







(16) 


where p^ Is the density of seawater. 


SOLUTIONS TO THE MODEL EQUATIONS 

, The time required to thin the lithosphere to a given fraction of Its 
Initial thickness was determined by Integrating (9) - (11) using a finite 
difference scheme. F1gu»*«« ^ gives relative thinning * fund on 

of the square-root of dimensionless time 

TS-^t , (17) 

P <=*0 

for two examples. Lithosphere thickness decreases with the square-root of 
time for a large fraction of the total thinning time. Near the equilibrium 
states, the /t - law Is not valid. (In the related classical Stefan problem 
of a freezing half-space the thickness of the solid layer Increases with the 
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squ&re>root of time.) We approximate as a function of time by 


0 i 


9S 


r — . o < t < (^) t^g 


*^0 “ ^f 


• ^ ^ ^95 * 


where t^g Is the time to accomplish 95X of the total thinning. Its dimension' 
less equivalent Tgg Is shown In Figure 2 as a function of the dimensionless 
parameters 

• <»> 


and 




( 20 ) 


The paran«ter X gives the total amount of thinning attainable and depends on 
three temperature gradients: the Initial subcrustal temperature gradient q^/k^ 
the final subcrustal t'imperature gradient q/k, and the slope. of the transition 
temperature profile. The parameter A is a measure of the total energy required 
for thinning fromt^j to It depends on the average temperature rise Tp(t*) - 
Tj^(t*) across the stagnation point boundary layer, on the transition temperature 
at surface pressure, on the temperature rise due to radioactive heat production 
In the crust, and on the ratio of the latent heat to the heat capacity. Figure 2 
shows that the thinning time Is proportional to A and to (1 - X)”, with n slightly 
larger than 2. 

Figure 3 also gives the ratio q* of the increase in surface heat flw 
q^ - (q^ is surface heat flow and q^^ Is the Initial value of q^) to 

“ ^s1 increase in surface heat flow after equilibrium has been re- 
established. The ratio q* increases from 0 to 1 as the lithosphere thins from 

^ • 
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Its Inltiil to Its final equilibrium thickness. The change In surface heat 
flow lags behind thinning by an amount which Increases with q/q^. At the time 
the surface heat flow has Increased by 10% of q^. • q^^ the lithosphere has 
thinned by 50% for q « 2q^ and by 90% for q * lOq^. 

^ % 

DISCUSSION 

We now apply our model to thinning of the continental and oceanic 
lithospheres on Earth and to thinning of the lithospheres of Venus and Men,. 
Geological observables will be calculated and inferences from the model cal- 
culations will be discussed. 

It is assumed that the material parameters pj^, o, c» k, and L take 
the sanrn values for each planet and tectonic region. These values are 
^to “ ® c - 750 J kg'^ k « 2.5 W ro’^ 

and L = 4.2 X 10® J kc]“^ for cmnplete meriting. They have been adopted from 
Stacey (1977). 

Two transition temperature profiles will be considered. In the first, 
the transition temperature is the 1000®C isotherm (models a). This is often 
thought to mark the transition from quasi-rigid to viscous behavior of mantle 
rocks over geological time (Oxburgh and Turcotte, 1978; Schubert, 1979) and 
it is also a reasonable approximation to the wet solldi of mantle rocks from 
laboratory experiments (Spohn and Schubert, 1982b). In the second, the t» in- 
sition temperature profile (modeli; b) is 0.75 times a representative dry 
solidus for the mantle given by Chapman and Pollack (1977). Here T^ « 1100®C 
and » 3.1 K km”^. 

We assume that the plume temperature Tp Increases with depth at a rate 
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of 0.4 K kro'^ which Is approximately the adiabatic temperature gradient In 
the upper mantle. Ue vurther rssume that Tp(i^) Is 100 K above the transi- 
tion temperature The average over of then takes 

values between 1(X} K and hOO K. 

Finally, we must consider what are reasonable values for tiie increased 
basal heat flux q. Heat ilow data for both continents and oceans on Earth 
vary by roughly a factor of 10 (Sclater et al., 1980). As Is evident from 
(4) heat flrw from the asthenosphere contr1bi«tes only partly to the surface 
heat flow. A factor of ten Increase In basal heat flow changes the surface 
heat flow by only a factor of five If Is approximately equal to DA^. Hantle 
plumes can transfer heat at significant rates. A plumj with flow velocity of 
10 mm yr*^ and tenperature 100 K hotter than the base of the llthospnerc ad- 
vects heat at the rate of about 1 W m*^, about 50 times the basal heat flow In 


low heat flow provinces. As discussed In the Introduction, this high heat 
flux might never be observed at the surface In the form of a conductive heat 


flow. The pardiT»ters q^ and depend on the planet and the tectonic region 


under consideration. For the oceanic crust Tj^ may be 25 K; for the continental 


crust It may be approximately 100 K. 


THINNING OF THE CONTINENTAL LITHOSPHERE 

Data on continental and oceanic heat flow have recently been reviewed by 
Srlater et al. (1980, 1981). Representative values of the continental heat 
flow arc close to 50 mW m' in heat flow provinces older than 250 Ma, while In 
younger provinces the mean Is 77 irtJ m*^. The scatter of data Is considerable 
in young provinces; Sclater et al. (1981) give a standard deviation of 53 mW m . 
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HoMever, It Is much smaller In provinces older than 2S0 Ha. Here Oie stin^rd 
deviation Is only 15 mW m"^. Taking ■ 50 mW m”^, ■ 2.5 wH m"^, and 

D ■ 10 km In (4), we get 25 mW m“^ as a representative value for the basal 
heat flux under continents. Equation (4) follows from (3). but It Is also 
Independently empirically valid (Birch et a1., 1968). It Is used to calculate 
a reduced surface heat flow from surface heat flow and radioactivity data. 

Sclater et al. (1981) give 25 iiAl m“^ to 29 inW ra*^ for the reduced heat flow 

from continents. The reduced leat flow may be a somewhat crude estimate of 

the basal heat flow because there Is radioactive heat production In the mantle. 

In addition, the continental lithosphere is still growing as a consequence 
of the cooling of the Earth and as a result the heat flux across the lltte- 

sphere Is some percent larger than the heat flux frcmi the underlying mantle 

(Schubert et al., 1979). 

The temperature rise across the crust due to radioactivity Is calculated 
from (6) to be 95 K, assuming a crustal thickness of 30 km. These values of 
and T|^ give Initial lithosphere thicknesses of 91 km (model a) and 146 km 
(model b) depending on the transition temperature profile assumed. The smaller 
value of corresponds to the lOOO^C Isotherm. These values are on the low 
side of recent assessments of continental lithosphere thickness. It has been 
argued that the continental lithosphere might be more than 300 km thick under 
old shields (e.g. Jordan, 1975; Sipkin and Jordan, 1975; Chapman ana Pollack, 
1977). However, Sclater et al. (1980) cohcluded that the continental litho- 
sphere thickness should not significantly exceed 150 km, an estimate In accord 
with other studies (Sleep, 1971; Crough and Thcxnpson, 1976b; Turcotte and McAdoo, 
1979; Bld>‘kweli and Chockal Ingam, 1980). 

In order to determine the circumstances under which convective thinning 
01 she continental lithosphere can be effective, we have calculated the time 
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required to thin the lithosphere by 955, tg^> for lithospheres 50 to 300 In 
thick subject to heat flux Increases from Z^^ to 3Cq^. These results are 
given In Fig. 3. A few cases are reported In more detail In Table 1. Addi- 
tional examples can be found In Spohn and Schubert (1982b). We list ^95 
whenever the completely thinned lithosphere Is thicker than the crust. In 
cases where the final thickness Is less than the cnistal thickness we give the 
time to thin to the crust as obtained from (18). The time to thin to 50 km 
thickness Is also listed for those cases In which the final equillbrim thictoiess 
Is smaller. The values of A vary from 0.11 to 0.50. We have assumed that 
not more than Incipient melting Is Involved. 

Very thick lithospheres (t^^ * 300 km) can be thinned to the crust In 
tens of minion years, the time scale usually associated with rifting (Seldler 
and Jacoby, 1981), If basal heat flow Increases by more than a factor of 10. 
Thinner lithospheres require less of an increase In basal heat flux, a factor 
of five can suffice. Because the radioactive heat production In the crust 
does not change substantially with time on this time scale, the increase In 
surface heat flow Is less than the increase in basal heat flow. Surface heat 
flow data from the rift valleys of contemporary continental rifts have been 
compiled by P. Morgan (1981). The values range frwn 70 to 125 mW m . This 
Is somewhat less than the surface heat flow of 139 mW m which obtains upon 
thinning the 100 and 150 km lithospheres of Table 1 with a factor of 5 Increase 
In basal heat flow. In these cases the lithospheres are thinned to 50 km In 
only 19 and 37 Ma. The smaller observed heat flows can be understood by the 
time lag between the increase in surface heat flow and thinning (Fig. 1) and 
by the consumption of part of the advected heat by partial melting once the 
lithosphere is substantially thinned. This heat will eventually reach the sur- 
face via volcanism. 
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The uplifts given In Table 1 correspond to lithospheres thinned to 
their crusts or to their final equlllbrluai thicknesses, idilctiever are larger. 
Uplift increases with 1^; It also Increases with q/q^ provided that If Is 
larger than It ranges from 1 or 2 km for low q/q^ and small to 

several kni for t1« larger values of q/q^ and The uplifts given here 
differ from those given In Spohn and Schubert (1982b) because here we have 
considered Isostatic equilibration with as the depth of compensation while 
In the earlier paper we considered uplift solely due to thermal expansion. 

Seldler and Jacoby (1981) give topographic heights of rift shoulders for many 
conteir^rary rift structures. With the exceptions of the Godawarl rift, the 
Hon Graben, and the Rhine Graben, where uplifts are less than 1 km, uplifts 
range from 1 to 4 km. These values agree with the ones given In our earlier 
paper and those given here for < 150 km. We have not accounted for erosion, 
which can occur at rates comparable to the calculated average uplift rates. 

For example, with = 150 kia and q » 5q^ we get an average uplift rate of 
0.06 km Ma”^. Allowance for erosion would reduce the calculated uplifts and 
Improve the fit to observation. However, our results argue against a 300 km 
thick continental lithosphere on Earth. It would not be easy to explain why 
plixnes would avoid such thick parts of the lithosphere If they existed. However, 
Mars' lithosphere Is likely to be around 300 km thick (see below) and uplifts 
around 10 km art* observed at Tharsis and other volcanic regions. 


VULNERABILITY OF THE CONTINENTAL LITHOSPHERE 

Our model does not account for the lateral m'^vement of the lithospheric 
plate with respect to the hot spot bel(».. It has been argued by Gass et al. 
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<1978) and by Pollack tt al. (1981) that the heat carried to the plate Ml 9 ht 
be swept away by plate notion thus precluding significant thinning. Pollack 
et al. (1981) define a 'vulnerability* paraneter V 

V r j"* »o . (S 

* j 

where u Is the plate velocity and d Is the characteristic spacing of hot 
spots. Essentially, they compare the propagation length (dk/itupc)^ In time 
d/u of a thermal wave to the Initial thickness of the lithosphere. This 
vulnerability parameter Is derived from the conductive model of Gass et al. 
(1978) which follows the upward propagation of a temperature disturbance at 
the base of the lithosphere. While the vulnerability paraneter varies from les! 
than 1 to 14 for the Earth's lithosphere. Pollack et al. (1981) show that SOX 
of the Earth's hot spot population Is located where V Is less than 5, I.e. 
where the propanatlon length of the thermal wave in time d/u Is more than 20X 
of the Initial thickness. 

We argue that the lithosphere should be more vulnerable to convective 
thinning than to conductive thinning. The relevant convective time scale ^95 
can be two orders of magnitude smaller thmi the conductive time scale for 
X > 0.85 (Fl'j. 2). Thus the amount of thinning attainable In d/u should be 
considerably larger for convective models as ccmipared with cwiductive ones. 
Taking 10 m as an order of magnitude value for d and 10 m a as one for 
u we fin'i that d/u Is 10® years. This is a factor of 10 larger than the time 
scale we obtain for our successful models. Moreover, the lithosphere is 
thinned by more than 50% in this time The finding of Pollack et al. (1981) 
that hot spots are preferentially located In regions where V is relatively 
small together with the ten million year time scale associated with mid-plate 
volcanism and rifting suggest that the heat flux from asthenospheric hot spots 
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Is not highly onoonlous. This Is because inw «fou1d expect the hot spots to be 
■»re evenly distributed over the Earth's surface If plumes advected much more than 
ten times the basal heat flux. The variations of heat flux Into the lithosphere 
need not be much larger than the variations In tN surface heat flux which are 
generally less than a factor of 10 (Sclater et al., 1981). 

£ 

THINNING OF THE OCEANIC LITHOSPHERE 

Heat flow data from the continents suggest that the mean basal heat 
f1(w Is approximately constant for provinces older than 250 Ha. ftean oceanic 
surface heat flow values from wll sedln^nted areas » on the contrary, decrease 
monotonlcally with crustal age. The data compiled by Sclater et al. (1980, 1981) 
fit the empirical relation 

[raw ra'^] = 509.3 t’** [mb] (22) 

for crustal ages less than 200 Ha. If we assume that the upper 5 km of the 

oceanic lithosphere Is a basaltic crust with a surface radioactive heat 

production rate per unit volume of 0.5 pW m and a characteristic depth 

of heat production D equal to 5 km, we get a contribution to due to radlo- 
_2 

activity of 5 raw m . The basal heat flow q^ Is obtained by subtracting 
5 mW m from q^. The temperature increase Tj^ due to crustal radioactivity 
Is 25 K. 

Again, we use the lOOO^C Isotherm and 0.75 times the generalized dry 
solidus from Chapman and Pollack (1977) to denote the lithosphere base. Initial 
lithosphere thickness as a function of age and q^ Is given In Fig. 4. At 
age 50 Ma lithosphere thicknesses are 36 and 43 km for the respective transition 
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temperature models. The thicknesses Increase to 79 and 111 km at an age of 

200 Na. For transition temperature model a (lOOO^C Isotherm), the amount of 

lithospheric thinning depends only on q/qp I.e. X ■ (q - q^)/q. Thus, the 

lithosphere can be thinned by 9S% If q * 20q^. by 90X for q ■ lOq^, by 80X f or 

q ■ 5q^, and by 50% for q ■ 2q^. For transition temperature model b, k * 

-2 

7.8 mW m and X depends on both q/q^ and q^ (see (19)). The dimensionless 
amount of thinning is 0.91 for q = lOq^. and q^ = 67 mW ro“^, the basal heat flux 
Into 50 Ma old lithosphere. This Increases to 0.92. for q^ « 31 nM m"^, the 
basal heat flux Into 200 Ma old lithosphere. For q/q^ = 5, X Increases from 
0.82 at 50 Ma to 0.84 at 200 Ma, while for q/q. * 2. X increases from 0.53 to 0.57. 
If q = 20 q^, X is 0.96. The parameter A is roughly constant for model a 
at the value 0.10. For model b, A increases with both and q/q^. Its 
maximum value here is 0.23 for * 111 km and q/q- = 20 and its minimum value 
Is 0.13 for * 43 km and q/q. = 2. Again, we have assumed that not more 
than incipient melting is involved. 

Figure 5 gives tgg, the time to accomplish 95% of the total thinning, 
as a function of the Initial thickness and q/q. for both transition tempera- 
ture models. The thinning time increases with and decreases with q/q^. It 
Is essentially independent of the particular transition tanperature profile. 

A comparison with Fig. 3 shows that tgg is also essentially independent of 
the thickness of the crust and the temperature increase due to radioactivity 
across it. A 100 km thick oceanic lithosphere can be thinned by 90% of its 
Initial thickness in 16 Ma if the basal heat flux is Increased by a factor of 
10. This is the same as the time needed to thin an equally thick continental 
lithosphere subject to the same heat flux increase. 

Taking a different approach, one may study lithospheric thinning over 
a plume delivering a constant heat flux q to the base of the oceanic lithosphere. 


The ratio q/q^ Is then a function of the age of the lithosphere rather than 
a constant paranieter as In the previous discussion. Recalling that * (1 - 
with n 5 2 and (1 - X) » we find that 

Tgg - q* " constant (23) 

for variations of (1 - X) that are not too large. As an example we choose 
.0 

q * 250 mW m » which Is the heat flow Into 4 Ha old lithosphere. It also 
Is about 5 times the heat flow Into 100 Ha old oceanic lithosphere and ten times 
the heat flux Into the continental lithosphere as discussed In the previous 
section. For ages of the oceanic lithosphere between 50 and 200 Ha, q/q^ 
then varies between 3.7 and 8.1. The oceanic lithosphere overriding such a 
plume win be thinned by 95X of Its initial thickness In about 10 Ha. 

Oceanic plates generally move faster than do plates that are largely 
made up of continental lithosphere. Thus the time scale d/u may be up to an 
order of magnitude smaller for oceanic plates than for continental plates 
(Pollack et al., 1981). Thinning of the oceanic lithosphere must occur at 
lease as fast as d/u. In the case of the Hawaiian Island chain the plate 
velocity Is 100 mm a“^, and d/u Is about 10 Ha. Crough (1978) has shown that 
the Hawaiian swell did rise to Its maximum height In 6 to 12 Ha. Slower moving 
plates like the Atlantic plate may not require such fast thinning times. It 
should be noted that Hawaii Is anomalous because It Is located In a region 
less vulnerable to thinning than the locations of 96% of all hot spots 
(Pollack et al., 1981). We conclude that the oceanic lithosphere, like the 
continental lithosphere, can be thinned In reasonable times by plumes deliver- 
ing less than 10 times the undisturbed basal heat flux. In fact, summarizing 
the results for the continental and oceanic lithospheres we note that mantle 
plumes need not advect more than 250 mW m“ . This Is 25% of the heat carried 
by steady state models of mantle plumes (Parmcntler et al., 1975; Yuen and 
Schubert, 1976). 
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UPLIFT OF THE OCEANIC LITHOSPHERE OVER HOT SPOTS 

Uplift due to thennal expansion and Isostatic equilibration of the 
oceanic lithosphere and replacement of lithosphere by hot asthenosphere Is 
given In Fig. 6. It Is assumed that the llthosphere-asthenosphere boundary 

# 

has reached Its final equilibrium position. The depth of Isostatic compen- 
sation Is the Initial lithosphere thickness 1^. Uplift Increases with q/q^ 
and with t^. The uplifts range from about 500 m for q/q^ equal to 2 and 
young oceanic lithosphere to 3 km (iiK}del a} and 4 km {model b) for 200 Ma old 
lithosphere and q/q^ » 20. Uplifts for q/q^ » 10 and q/q^ ■ 20 differ by no 
more than 5X; they are approximately twice as large as uplifts for q/q^ » 2. 

The dotted line In Fig. 6a gives the uplift for the model that has a plume 
advecting heat at a constant rate of 250 mW m . This uplift Is close to that 
calculated for q/q^ * 5. It is smaller for ages less than 75 Ma and larger 
for older lithosphere. 

Also given in Fig. 6 is subsidence as a function of age of the oceanic 
lithosphere relative to the height of the mid-oceanic ridge. This subsidence 
Is a consequence of the cooling and thickening of the oceanic lithosphere as 
it moves away fr«n the spreading center. On calculating the subsidence curves 
we have assumed that at the spreading center = 0 km and that the astheno- 
sphere there is at Tp. From Fig. 6 it is seen that the reheating and thinning 
of the lithosphere over a hot spot fully compensates the subsidence if the 
increase in basal heat flux is at least as large as a factor of 10. If q/q^ ■ 5, 

less than 400 m would remain uncompensated. Uplift due to reheating and 

-2 

thinning over a plume advecting 250 mW m would leave about 300 m of net 
subsidence. 

Crough (1978) has compiled data on mid-plate swells. He points out that 
the crests of contemporary swells rise to a roughly uniform depth of 4250 m 
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below sealevel, about 1750 m below the crests of the mld-oceantc ridges. 

This can be explained by postulating different degrees of partial melting 
beneath the mid-oceanic ridge and the swells. Assume that the density beneath 
the mid-ocean ridge and the swell Is except for a layer of thickness- a 
beneath the ridge where the density Is less than Pp by an mnount Ap^ due to 
partial melting. If As Is the height of the ridge above the swell then 


AS « a 


Ip 



(24) 


In Fig. 7 we present the data compiled by trough (1978) together with model 
determinations of the depths to mid-plate swell crests as a function of 
crustal age. The model assumes that the crust at the mid-oceanic ridge Is 
underlain by a 35 km thick layer of asthenosphere in which the density Is 
' 3 to 65J lower than density at the same depth beneath the swell. This density 
difference can be due to a higher degree of partial melting In the source 
region of mid-oceanic ridge basalts (MORB) compared to the source region of 
oceanic island basalts (OIB). It is known that OIB are often richer 
In alkalis than MORB which can be due to a lesser degree of partial melting 
of the source rock. 


THINNING OF THE LITHOSPHERES OF VENUS AND MARS 

Bi'cause the thermal evolution of a planet is governed by its surface 
area to volume ratio, lithosphere thickness should be Inversely correlated 
with planetary size (Kaula, 1981). This suggests that Mars' lithosphere Is 
thicker than both Venus' and Earth's. Thermal history calculations of litho- 
spheric growth have yielded a thickness of about 300 km for Mars (Schubert et 


aU, 1979). The same thickness has been derived from Isostatic models of the 
Tharsis region (Sleep and Phillips, 1979). If the Venusian mantle has the same 
rheology as the Earth's mantle and If the crusts of the two planets are similar 
In thickness and radiogenic heat production, then Venus' lithosphere would be 
thinner than Earth's because of Venus' high surface t«iperature. ikn«ever. It 
Is also possible that just the opposite Is true. Comparative studies of 
Venus' gravity field and surface topography show that they are highly corre- 
lated (Phillips et al., 1981; Reasenberg et al., 1982} and that Venus' topog- 
raphy may be less compensated than Earth's. Thus, Venus' lithosphere might 
be thicker and possibly more rigid than Earth's despite the higher Venusian 
surface temperature. This could be the case if Venus' mantle Is very dry 
compared with Earth's either because Venus had less water than Earth to begin 
with or It lost its Initial water content early In its evolution as a consequence 
of a runaway greenhouse. 

Crustal thicknesses on both Mars and Venus could be larger than on Earth. 

The crusts of these two planets may be as thick as 100 to 150 km (Wood and 

Head, 1978; Frey, 1979; Anderson, 1979; Sleep and Phillips, 1979; Phillips 

and Lambeck, 1980; Phillips et al., 1981; Schaber, 1982). Surkov (1982) has 

summarized the radiogenic element concentration of Martian and Venusian surface 

rocks Inferred from In situ analysis and remote sensing. There are rocks In 

* 

the surfaces of both planets with radiogenic element concentrations similar to 
those of basalts and granites on Earth although the Martian and Venusian 
concentrations are son^what lower than typical Earth values. If the planets 
have the same abundance of heat sources per unit mass then radiogenic heat 
production scales with their mass. Consequently, If Venus and Mars are not 
significantly more differentiated than Earth, they should have lower abundances 
of crustal heat sources integrated over crustal depth than Earth, regardless 
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of their crustal thickness. The temperature rise across the crust T|^ of Mars 
and Venus due to radiogenic heat should then not exceed the value for the 
continental crust on Earth which Is ~1(M) K. Vte can account for Venus surface 
temperature of 475^C by adding this temperature to 

There Is no reason to require lithospheric thinning on Mars to have 
occurred on the relatively fast 10 Ma timescale that Is representative of 
rifting on Earth. There Is evidence In the cratering record that the Tharsis 
region may be 3 Ga old and that volcanic activity has continued to the present 
or recent past ('Vise et al., 1979). Other volcanic provinces like Elysium, 

Hellas, and Hesperia may also be 1 to 4 6a old (Carr, 1976). The time scale 
associated with the growth of these regions could thus be as long as 1 Ga. 

On the other hand, a shorter time scale Is certainly possible and Carr (1976) 
has suggested that Olympus Mons on Tharsis might have formed In only 200 Ma. 
Assuming » 300 km, « 30 km, Tj^ * 100 K and all other parameters equal to 
their corresponding terrestrial values, we find that Mars' lithosphere can be 
thinned In hundreds of Ma by relatively modest Increases In heat flux, I.e. 
q/q^ < 10 (Table 1, Fig. 3). If the Martian crust were 100 km thick, the 
lithosphere would be thinned to the crust in less time. Even smalT.*r increases 
In heat flux would be required if thinning occurred early in the thermal history of 
Mars. Thermal history calculations suggest that the Martian lithosphere 
3 Ga ago, for example, was only about 100 km thick (Schubert et al., 1979). 

I . 

It may have remained thinned under old volcanic provinces ever since their 
formation. 

Turning to Venus, we note that for given q. and q/q. thinning would be 
faster than on Earth if Venus' model parameters are the same as Earth's. The 
reason is that Venus' lithosphere thickness would be smaller In this case 
because of its high surface temperature. According to Schaber (1982), the 
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excellent preservation of Venusian rift valleys and the helpht of their rlns 
of about 2 km suggest that rift structures on Venus are young and form on a 
time scale comparable to that of continental rifting on Earth. If Venus has 
a thick lithosphere with > 100 km, large Increases In heat flux (q/q^ > 10) 
would be required to thin the lithosphere In tens of Ha. On the other hand, 
Phillips et al. (1981) note that the Equatorial Highlands within which the 
rifts are located may be old features. Thus the argimients put forward for 
Mars ma> also apply for Venus. The lithosphere under the highlands may have 
been thinned by plumes of modest strength over times of hundreds Ha or even Ga 
and it may have remained thinned to the present day. (Note that the decay 
time of the lithosphere anomaly once the plirnie 1s switched off Is much larger 
than the thinning time (Spohn and Schubert, 1982b)). Recent rift structures 
may then be caused by further thinning and disrupting the pre-thinned lithosphere. 

Uplifts for models of lithosphere thinned to the crust on Mars and 
Venus are given In Tables 2 and 3. The crustal thicknesses are also the final 
equilibrium thicknesses of the lithosphere. The tables Include uplifts for 
the crustal thicknesses of 30 and 100 km used above. Uplift Increases with 
Initial lithospheric thickness and with q/q,|. It also Increases with 
decreasing z^. Uplifts are generally much larger for transition temperature 
models b, partly because Is larger for a given and partly because 
Tp(t*) - To(t*) Is larger. (Remember that T^ Is chosen so that Tp(*^Q) * 

T^(to) * 100 K). The larger value of q^ Increases (13). The larger 
value of Tp(t*) - T^(l*) Increases the absolute value of App (14). In the 
models presented here APp does not exceed 3% of pp. 

Martian uplifts In Table 2 range from 1 to 11 im. They would be about 

“5 -1 

1.5 times larger If the thermal expansion coefficient of 4.5 x 10 K were 
used as In some of the terrestrial models. A further Increase of uplift Is 
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possible If the plum Is partially melted. Density differences due to partial 
melting can easily Increase the absolute value of by another 3X of Pp. 
Hence* the model may explain the height of Tharsis and other volcanic regions 
of Nars of 10 km above the planet's median level (Carr, 1976; Greeley and 
Spudls, 1981). The large Martian uplifts are mainly a consequence of Mars' 
thick lithosphere. 

For Venus we calculate uplifts ranging from SOO m to 4 Imi. These match 
*«el1 the topographic heights of the ridges on Venus that rise 500 m to 2.5 kn 
above the surr<Hind1ng plateaus (Schaber, 1982). Uplifts on Venus are lower 
than uplifts on Earth and Mars for given and q/q^. The high surface 
tw^rature which adds to Tg reduces both and the absolute value of Opp 
because q^ Is smaller (5) and, more Importantly for thin crusts, T^(z) lies 
closer to the plume temperature profile. . Uplifts on Venus would be as large 
as uplifts on Earth and Mars lf the 11thosphere*asthenosphere transition 
temperature in the Venusian mantle were about 500 K larger than In both the 
terrestrial and Martian mantles. 

We conclude from our calculations that the model of convective thinning 
of the lithosphere can also be successfully applied to structures on Mars and 
Venus that resemble hot spot regions on Earth. 
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TABLE 2. Uplift for models of the Martian lithosphere thinned to 

the crust. Models a use the 1000®C isotherm to denote the 
lithosphere base while models b use 0.75 times the dry 
solidus. 
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TABLE 3. Uplift for models of the Venusian lithosphere thinned to the 
crust. Models a use the 1000®C isotherm to denote the litho- 
sphere base while models b use 0.75 times the dry solidus. 

T|^ is 575°C, all other parameters are given terrestrial values. 
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FIGURE CAPTIONS 

Fig. 1 Relative thinning as a function of the square^root of dimensionless 
time T s (k/pct‘)*t for Increased basal heat fluxes q • 2q^ and 
q • lOq^ (frnn Spohn and Schubert, 1982b). The relative thinning of 
the lithosphere as a function of tine can be approximated by the In- 
dicated ramp functions (light lines). The dashed lines give q^ h 
^** s ~ " ^s1^ where q^ is the surface heat flow, q^^ Is the 

Initial surface heat flow, and q^^ Is the surface heat flow after 
equlllbriimi is re-established. Note that q* lags considerably behind 
the thinning. 

Fig. 2 Dimensionless time = (k/pcijj)'tgg to achieve 95S of thinning as 

a function of relative thinning X = (£.^ - (a) and t (b). 

Fig. 3 Thinning of the continental lithosphere. Time tgg to achieve 95% of 

thinning versus Initial thickness The curve parameter Is the basal 
heat flux ratio q/q^. Circles refer to transition model temperature a 
wherein 1000°C denotes the base of the lithosphere. Crosses refer to 
transition temperature model b in which the lithosphere base Is denoted 
by 0.75 times a representative dry mantle solidus. 

Fig. 4 Thickness of the oceanic lithosphere versus age and surface heat flux. 
Model a uses the lOOO^C isotherm to denote the lithosphere base while 
model b uses 0.75 times a representative dry mantle solidus. 

Fig. 5 Thinning of the oceanic lithosphere. Time tgg to achieve 95% of 

thinning versus initial lithosphere thickness with q/q^ as a paraiwter. 
Circles refer to transition model temperature a and crosses to model b. 

The triangles are from the continental lithosphere models of Fig. 3. 
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Fig. 6 Uplift upon thinning of the oceanic lithosphere for transition temperature 

models a and b as a function of age and Initial lithosphere thickness 

with heat flux ratio q/q^ as a parameter. Here models a and b differ 

also In the volumetric coefficient of thermal expansion which Is 4.5 x 10’^ 

for model a. Also given (dashed lines) Is the subsidence due to cooling 

of a column of undisturbed oceanic lithosphere as It moves away from the 

mid-oceanic ridge. With q/q| > 10 uplift upon thinning can fully compen- 

sat this subsidence. The dotted line refers to a model in which the 

-2 

plunm advects heat at the constant rate of 250 mW m . 

Fig. 7 Water-depth to the crests of mid-plate swells versus crustal age. The 

data are from Crough (1978). The crests rise to a common depth of about 

4 km. The thin lines indicate the results of model calculations that 

« 

explain this observation. The model assumes that the crust at the 
mid-oceanic ridge Is underlain by a 25 km thick layer of asthenosphere 
In which the density is 3 to 6 X lower then the density at the same 
depth beneath the swell. SH: St. Helena, K: Kerguelen, M: Marquesas, 

CA: Cook-Austral, R: Reunion, S: Society, H: Hawaii, T: Trinidade, 

B: Bermuda, CV: Cape Verde. 
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